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ABSTRACT 

We present a systematic study of stellar feedback processes in simulations of disk 
galaxy formation. Using a dark matter halo with properties similar to the ones for the 
Milky Way's stellar halo, we perform a comparison of different methods of distribut- 
ing energy related to feedback processes to the surrounding gas. A most promising 
standard model is applied to halos spanning a range of masses in order to compare 
the results to disk galaxy scaling relations. With few exceptions we find little or no 
angular momentum deficiency for our galaxies and a good agreement with the angular 
momentum-size relation. Our galaxies are in good agreement with the baryonic TuUy- 
Fisher relation and the slope of the photometric TuUy-Fisher relation is reproduced. 
We find a zero-point offset of 0.7 to 1 magnitudes, depending on the employed IMF. 
We also study our standard feedback model in combination with additional physi- 
cal processes like a UV background, kinetic feedback, a delayed energy deposition as 
expected for type la supernovae, mass return and metal-dependent cooling. Only a 
combination of effects yields a real improvement of the resulting galaxy by reducing 
the bulge, while including metal-dependent cooling increases the bulge again. We find 
that in general the stellar mass fraction of our galaxies is too high. In an ad-hoc ex- 
periment we show that an removal of the bulge could reconcile this. However, the fit of 
the TuUy-Fisher relation can only be improved by delaying the star formation, but not 
suppressing it completely. Our models do not seem to be efficient enough to achieve 
either effect. We conclude that disk formation is a complex, highly interconnected 
problem and we expect a solution to come from a combination of small effects. 

Key words: galaxies: spiral - formation - evolution - structure - methods: A^-body 
simulations - hydrodynamics 



1 INTRODUCTION 

The formation of galaxies is a complex process governed 
by gravitational collapse and an energetic coupling between 
the interstellar (and also intergalactic) gas and stellar evo- 
lution processes. This includes the cooling of gas and sub- 
sequent formation of stars as well as loss of mass and en- 
ergy through stellar winds and supernova explosions which is 
known as stellar feedback. Such feedback is essential for the 
success of the widely accepted hierarchical model of galaxy 
formation (fwhitc & Rees 1978). It regulates star forma- 
tion by particularly preventing stars from forming too early 
which would lead mostly to g alaxies dominated by large, 
old spheroids instead of disks (|White fc FrenkHlOQll ). This 
is inconsistent with observations showing that up to 70% of 
10^^ h~^ Mq halos host disk do minated late-type galax- 



gas and disrupt cold clouds, therefore efficiently quenching 
star formation. At the same time, the surrounding gas is 
enriched with metals produced in stars ena bling more ef- 
ficient cooling l|Sutherland fc Dopita I Il993t) and therefore 
increasing star formation. All of these processes happen 
on scales of a few to tens of parsecs which is well be- 
low current achiev able resolution of c osmological simula- 
tions (see however ICeverino fc Klvpini ,2009 ) . Therefore a 
range of so called "subgrid models" have been developed 
to model at least the ef f ects o n larger, ga lactic sca l es (fo r 
example lOkamoto et ahl (|2005|), [Sc annapi eco et al. I l|2006l ). 
ISpringel fc Hernguistl (|2003l ). IStinson et al.1 ||2006|)'). Com- 



ies in the presen t day universe (e.g. IWeinmann et alTl2006l : 
iPark et al.ll2007l ). Feedback affects the thermodynamics of 
the interstellar medium (ISM) in two ways. In supernova 
explosions, large amounts of energy heat the surrounding 



pared to early work (for exa mple ISteinmetz fc Miilled 
1 19951 : In avarro fc Steinmetzlll997h and in combination with 
improved resolution, this has recently led to consid- 
erable improveme nt in reprod u cing individual realistic 



is 



Abadi et al 



disk g al axies , . .. .. . _ 

||2004). iGovernato et al.1 ||2007|'). lOkamoto et a.1 



Governato et al.l 



[Robertson et al.1 (|2004l '). IScannapieco et al.l (|2008l ')'). though 
still none of them fulfilled all characteristics of typical ob- 
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served late-type spirals. Furthermore, none of them resem- 
bled a bulge-less disk galaxy, de spite of their non-neg ligible 
numbers in the local universe ()Kautsch et al.ll2006l ). The 
fundamental question of disk galaxy formation therefore is 
still open. 

All successful attempts to model disk galaxies which 
have been achieved so far were simulated using handpicked 
dark m atter halos and ca refully calibrated codes. Recent 
work bv lScannapieco et al.l (2009) demonstrated that, when 
applying their model to a random sample of halos, only a 
small fraction had a disk component and n o ne w as disk 
dominated. Furthermore, IScannapieco et ahl l|2009l ) ques- 
tions the usual assumption that Hubble type and forma- 
tion history are directly correlated. Considering the large 
abundances of disk galaxies, this indicates that we may 
still lack an understanding of the big picture. Furthermore, 
disks are susceptible to mergers which are more common 
at high rcdshift and typica lly result in t he format ion of a 



spheroid (Peebles 1982: Blu men thal et al. [ l984; Davis et al 
' Gottlober et al.ll200ll : IWechsler et al.ll2002l : ICole et al 



1985 




2008 





Wetzel et al.ll2009l ). One possible way out is to re-grow 



a disk after a major merger, provided the merging progeni- 
tor d isks are gas-rich l|Robertson et a l. 2006; Bullock et al.l 
l2009i V A first example of this in a cosmological simulation 
has been presented recently bv lGovernato et al.l (|2008h . but 
again it is an isolated case. 



The work we present in this paper focuses on furthering 
the understanding of the physical requirements of disk for- 
mation. The main question addressed by our feedback study 
can be phrased as follows: Which of the many ingredients of 
the range of models summarized above are necessary, which 
are sufficient and how exactly do they influence the galaxy 
formation process alone and in combinati on? As in our study 
of th e angular momentum problem (Piontek fc Steinmet j 
I2OO9I . Paper 1 hereafter), we follow our strategy of slowly 
increasing the complexity of the model. First using a prese- 
lected dark matter halo with a mass similar to the halo of the 
Milky Way, with a high spin parameter and a quiet merging 
history, usually thought to be the most likely to host a disk, 
we study different methods of distributing feedback energy. 
The best result becomes our standard model, which we then 
apply to a set of halos with a range of different masses and a 
second set of Milky Way-type halos but with different merg- 
ing histories. We also study a set of simulations examining 
other relevant physical effects. We carefully analyze the re- 
sulting disk to understand the influence each of these effects 
has alone, before we put several physical effects together for 
our most realistic model of galaxy formation. 

The paper is organized as follows. In §2 we describe 
the initial conditions, the code and the analysis methods 
we use. The study of feedback energy distribution methods 
is covered in §3. §4 contains the study of halos with differ- 
ent assembly histories, §5 the study of halos with different 
masses and §6 the application of additional physics besides 
the standard feedback model. A detailed description of our 
most complex model can be found in §7 and the influence of 
the additional physics on the scaling relations is discussed 
in §8. We discuss our findings and conclude in §9. 



2 INITIAL CONDITIONS, CODE AND 
ANALYSIS METHODS 

2.1 Code 

The simulation s are done with the N-body code GADGET2 
(|Springelll2005l ) using t he Smoothed Particle Hydro d ynam- 
ics (S PH) framework (|Gingold fc MonaghanI (119771 ), ILucvI 
(|1977) )) for the gas and an implementation of ra diative cool- 
ing co urtesy of Volker Springel and based on iKatz et al.l 
(Il99d ). For simplicity, in the standard cases we do not in- 
clude an external UV b ackground. Star formation is imple- 
mented following iKat j (1992) and is already des cribed in 
detail in Paper 1. It is motivated by a Schmidt law (jSchmidtl 
I1959I ) giving the star formation rate as 



' U 



(1) 



with c^,=0.1 and t*=max(tdyn,tcooi), and then the applica- 
tion of a stochastic approach. We assume that each gas par- 
ticles can spawn two generations of stars, so each star parti- 
cle has half the mass of a gas particle. Necessary conditions 
for star formation are a critical density of pcrit=7x lO^^^'g 
cm""^, a converging flow and a low temperature of T< 3 x 10* 
K. 



2.2 Initial conditions 



All o f our halos are resimulations (see iNavarro" fc Whitd 
Il994h from a large cosmological box wi th 64 h~'^ Mpc^ box 
size. We use the WMAP3 cosmology (|Spergel et al.l |2007| ) 
with Ho=73 kms~^Mpc"\ cr8=0.75, nrms=0.9, ^0=0. 24 
and f2A=0.76. The simulations start at z=50. Gas particles 
are added on top of each dark matter particle assuming the 
cosmological baryon density, nbar=0.04. Our standard halo 
(MW_mr) is a halo with a mass s i milar to the Milky Wa y 
(about 10^^ h-i Mq ISmith et all l2007l : IXue et all l2008l l. 
with a fairly quiet merging history and a comparatively high 
spin parameter A=0.054. Our standard resolution corre- 
sponds to 1024^ effective particles within the high resolution 
region surrounding our target halo. This translates to parti- 
cle masses of 2.8 x 10® h"^ Mq for gas and 1.34 x 10^ h"^ Mq 
for dark matter particles, and we use a softening of 1.5 and 
2 h^^ kpc respectively. Additionally, we selected three halos 
comparable in mass to the standard halo, but with very dif- 
ferent merging histories, to study the influence of the merg- 
ing history of the host halo on the formation of the disk. 
Lastly, we selected six halos with smaller masses starting at 
~ 10^^ M© and one halo with a larger mass of « 2 x lO^'^ M© 
to sample the disk galaxy scaling relations over a larger mass 
range. All halos were selected to have no object of equal 
or larger mass close by. All the lower mass halos and also 
the standard halo have been simulated in high resolution 
with 2048^ effective particles in the resimulation region and 
particle masses and gravitational softening parameters of 
3.54 x lO'^h"^ Mq, 1.7 x lO^'h"^ Mq and 0.75 and 1 h"^ 
kpc, respectively. Table[T]shows an overview of characteristic 
parameters of the halos we simulate. 

2.3 Analysis 

We define our halos using Rvir via p{Rvir) = Apc with 
Pc = 3H^/(8tvG) and A = IStt^ + 82a: - 39a;^ x = 
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halo 


Mvir 
1012 Mq 


Kvir 

kpc 


Vvir 

km 


A 


(1:3) 


effective 
resolution 


iNvir 


MW_mr 


1.26 


296 


146.17 


0.054 


3.63 


1024^ 


1.6 X 10^ 


MW.hr 


1.23 


290 


143.22 


0.058 


3.63 


2048^ 


10® 


MW_mhl_mr 


1.24 


294 


145.66 


0.023 


4.85 


1024^ 


1.55 X 10^ 


MW.mh2_mr 


1.23 


298 


137.2 


0.063 


1.18 


1024^ 


1.8 X 10"' 


MW.mh3.mr 


1.2 


291 


144.2 


0.05 


1.35 


1024^ 


1.5 X 10"' 


DMJirl 


0.135 


139 


68.32 


0.04 


4.56 


2048^ 


10^ 


DMJir2 


0.252 


173 


85.11 


0.029 


3.02 


2048^ 


2.42 X 10^ 


DMJirS 


0.358 


193 


95.34 


0.034 


3.44 


2048^ 


3.13 X 10"' 


DM_hr4 


0.493 


214 


105.5 


0.019 


6.09 


2048^ 


4.2 X 10^ 


DMJirS 


0.594 


227 


112.38 


0.016 


3.44 


2048^ 


5 X 10^ 


DM_hr6 


0.703 


241 


119.36 


0.026 


3.83 


2048^ 


6.1 X 10^ 


DM_mr7 


2.36 


364 


180.37 


0.018 


4.556 


1024^ 


2.85 X 10^ 



Table 1. Characteristic parameters for our halos run with the standard feedback model. MW labels our standard, Milky Way- type 
halo, MW_mh halos with the same mass but varying merging histories (mh) and DM_hrl to 7 halos with different masses. Ir, mr and hr 
indicate the effective resolution of the high resolution region shown in the second to last column. Nvir is the approximate total number 
of gas, star and dark matter particles in Rvir in the runs with the standard feedback model. This slightly varies for different feedback 
models, but the order of magnitude stays the same. 



no(l + 2)^/(^0(1 + zf^ + Ha) — 1. Characteristic param- 
eters for the galaxy are calculated using all particles within 
a sphere of 6Rd. For the Milky Way with Rd ~3.5 kpc, this 
corresponds to the disk radius. For our simulated galaxies 
with larger scale lengths, 6Rd corresp onds approximate ly to 
the size of the Milky Way stellar halo ISchneideil l|200d) . Rd 
is the disk scaling length. As gaseous disk we define the cold 
gas (with r < 3 X lO^K) in the center of the halo, which 
forms a clearly distinguishable disk. In most cases, there is 
only little hot gas surrounding the cold disk. We rotate the 
galaxy so that the disk is in the x-y plane. Rotational veloc- 
ities are measured at 2.2Rd using the rotation curve of the 
disk gas. This is close to observational measures of rotational 
velocities wh ich are often obtain ed using HI observations 
(for example Ide Blok et ahllioOSl ). An example for rotation 
curves with the standard halo and standard feedback model 
is shown in Figure [T] for the standard and the high resolu- 
tion run. Our rotation curves are fairly fiat, though particu- 
larly in the smaller halos, and in the more bulge dominated 
models, they can be steeper than shown here. However, the 
difference between the measured rotational velocity at the 
peak and at the adopted radius, even for our steeper curves, 
is never more than 10% and therefore does not change our 
results significantly. Particularly it cannot account for the 
offset of the zero-point of the TuUy-Fisher relation we see 
for our simulated galaxies. 

The angular momentum of the galaxy is calculated for 
the gas and stars of the galaxy as described in Paper 1, 
by summing over all particles, with respect to the cen- 
ter of mass. This angular momentum is compared to the 
value expected from the properties of dark matter halos 
.icaic ^ 1-3 X 10^(Vro t /200)2 kms~^h'ikpc as derived by 
iNavarro fc Steinmet j l|2000bl ). 



In order to study the importance of the disk in our 
galaxie s, we also p e rform a dynamical decomposition fol- 
lowing lAbadi et al.) (|2003bh . We compare the z component 
of the angular momentum jz with the angular momentum 
of the corresponding circular orbit jcirc(E) for each star par- 
ticle. The ratio ej = jz/jcirc(E) describes the degree of rota- 
tional support of a given stellar particle. A thin, rotationally 
supported disk has ej ~ 1, a spheroid has little net rotation 
due to equal numbers of stars on co- and counterrotating 
orbits, its distribution of ej peaks at zero. All stars fitting 
in neither category are classified as thick disk stars. They 
are not rotationally supported at the same level as the thin 
disk, but rotate in the same manner, and spatially form a 
thick disk. The thin disk is generally dominated by young 
stars, but can also contain older stars. The dynamical bulge- 
to-disk ratios quoted in the tables are ratios of the mass of 
the bulge to the combined mass of thin and thick disk. 



2.4 Comparison with observations 

In order to perform a realistic comparison with observa- 
tional results, particularly with respect to observed scaling 
relations like the TuUy-Fisher relation, we create and an- 
alyze mock observation s following a method described by 
iKhalatvan et al."1 (|2008|). Usi ng population synthesis models 
bv lBruzual fc CharlotI (|2003| ) (BC) and the software SKY- 
MAKE^ we mimic observations performed by the Hubble 
Space Telescope. For this, the metallicity of each star par- 
ticle is needed. For simplicity, in our runs without explicit 



1 Berlin & Fouque (2007), 

http: //terapix.iap.fr/rubrique.php?id_rubrique=221 



4 F. Piontek & M. Steinmetz 



250 






1 1 1 1 1 


1 1 1 1 






1 1 1 1 1 1 1 1 1 








-- 










/ 














200 


- 1 
I : 






— 


— 








- 1 /f 
" f / 














150 


-llf 
""'/' 
"'/' 














100 


7/ 














50 






























potential — — — 

young stors 

disk qos 

2.2*R, 













,1 









5 10 15 20 25 30 35 

X [kpc] 

Figure 1. Rotation curves for tlie standard feedbaclc model in 
standard (red) and high resolution (black). The long dashed line 
dubbed 'potential' is the circular velocity curve calculated from 
Vrot = A^/ GM{< R)/R, the solid line traces the young stars in the 
rotationally supported disk, the dash-dotted line the cold disk gas. 
The vertical lines mark the location 2.2R(j, where the rotationally 
velocity is measured. 



metal enrichment, we assume solar metallicity for star par- 
ticles younger than 10 Gyr, and Z=10~ for star particles 
older than 10 Gyr. We explor e BC models w ith Salpeter 
(|Salpeteilll955l ) and Chabrier IChabried 12001 ) initial mass 
functions (IMF). The main difference between the two is 
a flatter (more physical) behavior at the low mass end for 
the Chabrier IMF, resulting in mass-to-light ratios being a 
factor of 1.5 smaller than for a Salpeter IMF. The images 
shown in this paper combine U, V and B band "data". We 
also create I band FITS files centered on the galaxy which 
are used to compute the surface brightness profile using the 
ESO-MIDAS pa ckagfl which is then fitted with GALFIT 
l|Peng et al.ll200^ . To obtain the galaxy's total magnitude as 
well as a photometric decomposition we fit a two component 
bulge-|-disk model to the whole profile. To find the disk scale 
length, we fit a pure exponential model to the outer parts. 
The rotational velocity is then measured with the rotation 
curve for the cold gas at 2.2Ri. The I ban d luminosity Li 
is calc ulated from the magnitude following ICourteau et al.l 
120071) as Li=10-°-^(^' 

Observers measure the angular momentum of a galaxy 
from the rotation curve and the scale radius as jobs = 
2_RdKot(2.2_Rd)- This relation, also known as the angular 
momentum-size relati on, is inferred for an exponential disk 
in an isothermal halo (|Fall fc Efstathiou 1980). For our sim- 
ulated galaxies, besides the real angular momentum content 
we also compute the angular momentum in this manner, 
using the I-band scale length and rotational velocity as de- 
scribed for the TuUy-Fisher relation. This "exponential disk 
estimator" is then used to compare to the actual angular 
momentum to investigate possible systematics. 

2 www.eso.org/projects/esomidas 



3 SUPERNOVA ENERGY DISTRIBUTION 
METHODS 

In cosmological simulations, a star particle basically corre- 
sponds to a single stellar population (SSP). The energy from 
supernova explosions is therefore summed over the whole 
SSP and then distributed over the surrounding gas particles. 
As a first step in our investigation of supernova feedback, we 
study three different ways of distribution. 

3.1 Implementation 

We assume a Miller-Scalo initial mass functi on for the stel- 
lar po pulation in one of our star particles l|Miller fc Scald 
(Il979h ) with a lower cutoff of 0.1 M© and an upper cutoff 
of 100 Mq: 

O.1<M<1M0 
l<Af<2M0 
20.3^-2.3 2<M<1OM0 (2) 

10 2°-^M-^-^ 1O<M<1OOM0 

where A=0. 284350751. Stars with masses between 8 and 40 
Mq explode as type II supernova, each explosion yielding an 
energy of 10*^ ergs. This results in an energy of 1.21x10*® 
erg per solar mass formed. This energy is smoothed over the 
neighboring gas particles using the SPH smoothing kernel, 
so each neighbor gets an energy of 

l^-HSN.i — CjSN,tot (O) 

P* 

/i* is the smoothing length and regulates the number of gas 
particles receiving feedback energy. For this we try two dif- 
ferent approaches. Initially, we use a fixed radius for the 
smoothing sphere of 1.37 h~^ kpc, which ensures a distribu- 
tion of energy consistent with resolution. However, in this 
scheme, the impact of the feedback energy is highly depen- 
dent on the local density, with the same amount of energy 
distributed to large or small numbers of particles. The struc- 
ture of the disks is improved when the size of the smooth- 
ing region scales with the local density instead, as does the 
smoothing length for SPH calculations in Gadget. With this 
scheme, the mass-weighted number of neighbors receiving 
feedback energy stays constant. 

In the following, we briefly describe the three meth- 
ods of supernova energy distribution. The instantaneous ap- 
proach is the simplest and also most drastic method. The 
total supernova energy is distributed to the neighboring gas 
particles immediately after formation of the star particle. 

Dumping the full amount of supernova energy in a single 
timestep is a very crude way of mimicking a real supernova 
remnant which expands and interacts with the surrounding 
ISM over an extended period of time, typically 20-30 Myr. 
The following two methods therefore distribute the energy 
more smoothly over time. In the heat-rate approach we dis- 
tribute the energy in equal portions over a total time of 
tsN=20 Myr by 

AE^EsN^, (4) 

tSN 

where At is the current timestep of the star particle. A sec- 
ond possibility for a slow energy output is the exponential 
approach, an exponential feedback model using 
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6.09 


224.8 


5.65 


2.41 


0.91 


0.92 


0.55 


MW_mrJ 


7.49 


8.1 


0.126 


0.173 


0.89 


7.9 


207.71 


2.63 


2.36 


0.31 


2.61 


0.95 


MW_mr_inc 


6.17 


7.53 


0.141 


0.127 


2.52 


5.85 


187.46 


3.27 


1.62 


0.57 


1.84 


0.9 


MWjnrJi 


10.7 


11.0 


0.136 


0.119 


1.12 


6.18 


239.37 


3.82 


2.52 


0.75 


5.06 


0.48 


MW_mr_hnc 


7.07 


8.91 


0.146 


0.079 


3.86 


5.36 


186.67 


4.23 


1.48 


0.42 


0.14 


0.86 


MW_mr_c 


10.05 


10.55 


0.134 


0.119 


1.51 


7.81 


221.06 


3.66 


2.47 


0.5 


2.48 


0.71 


MW_mr_enc 


8.05 


8.7 


0.136 


0.094 


1.31 


7.03 


221.2 


3.3 


2.21 


0.35 


0.59 


0.7 


MW_mr_encv 


10.4 


10.8 


0.143 


0.104 


2.55 


6.01 


243.1 


4.63 


2.08 


0.43 


0.27 


0.62 


MW_hr_cncv 


7.2 


7.82 


0.113 


0.11 


1.38 


7.89 


207.47 


3.13 


2.06 


1.02 


1.49 


0.9 



Table 2. Characteristic parameters for the galaxies with different distribution methods for the feedback energy. Vrot is measured for 
young stars at 2.2Rd,i. The luminosity and mass-to-light ratio are based on a Chabrier IMF for the edge-on projection. The last column 
measures the angular momentum content and contains the ratio of the total angular m omentum of stars and gas in the galaxy and the 
expected value based on the relation jcalc ~ 1-3 X 10'^(K-ot/200)'^ kms~^ h~^ kpc from lNavarro fc Steinmetj l|2000bl V 
^ galaxy mass including cold gas only, in 10^" Mq; galaxy mass including all gas, in 10^*^ Mq; ^ in WIq yr~^; * in kpc; ^ in km s~^; 
in IQlO L© 
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Figure 2. Mock observations in face on and edge on projections 
of our simulated galaxies with different methods of distributing 
feedback energy as described in the text and in Table |2] Each box 
has a side length of 60 kpc, and those with the edge-on projection 
a height of 30 kpc. In the top row we have the model without 
feedback, the instantaneous model, the heat-rate model and the 
exponential model. The models in the bottom row correspond to 
those on the top, but with cooling locally turned off. The bottom 
left panel is our standard model (exponential with local cooling 
turnoff and variable smoothing length). 



AE = EsN 



(5) 



t-t^ -inti At 
e *SN . 

tsN tsN 
Again, tsN=20 Myr and t* is the formation time of the stel- 
lar particle. This model mimics effects owing to the forma- 
tion time and lifetimes of stars with different masses. 

We also test variations of all of these methods where ra- 
diative cooling is turned off temporarily in the gas particles 



receiving supernova energy (iGerritsenI Il997l ) . This helps to 
prevent immediate re-radiation of the thermal energy due to 
the high density of the gas surrounding the star formation 
sites and the inability to resolve the multiphase structure of 
the interstellar medium. We choose a turnoff time period of 
20 Myr, based on typical cooling times for supernova blast- 
waves. With this we attempt to suppress the strong early 
star formation more efficiently. 

3.2 Results 

Table [5] and Figure [2] give an impression of the results for 
the different energy distribution methods. Our standard 
halo MW_mr has been used as initial conditions for all of 
these runs. The abbreviations i (^instantaneous), h (=heat- 
rate) and e (=exponential) indicate the type of feedback, nc 
stands for the local cooling turnoff, and nf for no feedback 
which is a run with only star formation which is included 
for comparison. All runs are employing fixed smoothing 
length except one, MW_mr_encv, with variable smoothing 
length for feedback energy. Variable in this case means that 
the smoothing length is calculated following the smoothing 
length for the SPH density calculation. It is scaled accord- 
ing to the mass- weighted number of neighbors which is kept 
constant to a value of 40. We experimented with the vari- 
able smoothing length also in the other cases but could not 
see substantial improvement. The first obvious result from 
Figure [2] is the importance of turning off the cooling locally 
in order to form a young, bright disk (bottom row vs top 
row). This is also clear from the stellar surface density pro- 
files shown in Figure [S] While all runs with feedback are 
quite similar outside of about 8 kpc, the galaxies with the 
cooling turn-off have a much flatter slope in the inner re- 
gions, corresponding to a smaller bulge. In the simulation 
without any feedback, the bulge is most dominant, as is ex- 
pected. The two prominent spiral arms seen in the image in 
the top left panel of Figure [2] are due to ongoing mergers of 
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Figure 5. Dynamical decomposition of the stellar component for 
the heat rate model (left panel) , the heat-rate model with variable 
smoothing length (middle panel) and the exponential model with 
exponential smoothing length (right panel), all with local cooling 
turnoff. The black line represents all stars, the blue line thin disk 
stars, the green line thick disk stars and the red line the spheroidal 
component. 



Figure 3. Stellar surface density profiles for the different feed- 
back models. The different colors indicate the different models 
as shown in the legend. The dashed lines show the corresponding 
model with local cooling turned off. The bump in the model with- 
out feedback and in the exponential model at about 25 kpc radius 
corresponds to an incoming satellite galaxy also clearly visible in 
Figure [2] 




Figure 4. Comparing the star formation histories for the different 
feedback energy distribution methods in comparison to the run 
without feedback (black line). The red line shows the standard 
case, the blue line the case with locally turned off cooling ("nc"). 




Figure 6. Left panel: Agreement of the galaxy formed in the 
standard high resolution model with the Kennicutt law, locally 
given by equation \E\ (solid line). In the outer region, the star 
formation rate drops since the disk is no longer Toomre instable. 
Right panel: Agreement of all our galaxies with different feedback 
distribution methods with the observed global Kennicutt l a w. Th e 
triangles are observational data points from iKennicut 
Our model galaxies populate the lower surface density region in 
agreement with the regular galaxies, while the observed galaxies 
with high surface densities are starburst galaxies. 



smaller satellites (one is also visible in the image) inducing 
new star formation. Nevertheless the galaxy is very massive, 
slowly rotating and clearly spheroidal with a very low angu- 
lar momentum. The clearly visible incoming satellite can be 
seen as a bump in the surface density profile. The bulge is 
also important in runs with feedback but without the local 
cooling turn-off. 

Figure |4] shows the star formation histories for all cases 
in comparison with the no-feedback run. We see that the 
instantaneous model even without the local cooling turn- 
off suppresses early star formation most efficiently due to 
the instantaneous energy input. The star formation history 
for the two models with smooth energy input are very sim- 
ilar, except when cooling is turned off locally, which has 
less influence in the exponential model than in the heat- 
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rate model. However, in the basic instantaneous model, star 
formation is suppressed but not self-regulated which is why 
we do not see a young stellar disk. The feedback actually 
becomes too strong. When we turn off cooling locally, a sig- 
nificant amount of gas stays hot for a longer amount of time, 
but also close to the disk, so when it cools again it can form 
stars there and we find a young but rather thick and unstruc- 
tured stellar disk. In the exponential model, the suppression 
is less efficient and the young disk less prominent. 

A more prominent spiral structure is achieved when a 
variable smoothing length is used, since in that case the 
code can better adjust to local density structures. For the 
exponential model with cooling turnoff this is shown in the 
bottom left image of Figure[2] Our two best cases in compar- 
ison to the Milky Way are MW_mr_hnc and MW_mr_encv, 
i.e. the heat-rate and the exponential model with local cool- 
ing turnoff and, for the latter, with a variable smoothing 
length. MW_mr_hnc has an almost perfect exponential sur- 
face brightness profile and is almost not angular momentum 
deficient. Nevertheless, as shown in Figure O the dynamic 
decomposition shows it to be rather a thick disk instead of 
a thin disk and a bulge, and the structure of the disk is 
very clumpy. This does not improve much when a variable 
smoothing length is used in this case. We therefore select 
MW_mr_encv as our standard model for the rest of our work. 
It fits the Kennicutt relation, given by 



EsFR = (1.5 ±0.7) X 10" 



M 



0pc- 



Mr. 



yr kpc' 



■,(6) 



well as shown in the left panel of Figure [HI though the outer 
regions are no longer Toomre instable and the star formation 
density drops. In general, all of our models fall well within 
the scatter of the observed global Kennicutt relation in the 
region of regular galaxies, as shown in the right panel of 
Figure [S] 

We test the standard model with our higher resolu- 
tion halo, as run MW_hr_encv, with about 10^ gas, star 
and dark matter particles within the virial radius. Our feed- 
back method does depend somewhat on resolution with the 
main effect being a more efficient suppression of star for- 
mation which reduces the total galaxy mass as measured in 
a radius of 6Rd by about 30%. The high resolution galaxy 
has basically no angular momentum deficiency. Even though 
the scale length in the high resolution case is a little larger, 
the surface brightness profiles are very similar. Overall the 
higher resolution is slightly different in some aspects but 
the model is reasonably robust and the changes at higher 
resolution improve the disk rather than having a negative 
impact. 



• 




(b) 

Figure 7. Panel (a): Mock observations of the Milky Way type 
galaxies with different merging histories, all run with the standard 
feedback model. For comparison the standard halo is shown on the 
left, followed from left to right by MW_mhl_mr, MW_mh2_mr and 
MW_mh3_mr. There is no edge-on projection for MW_mh2_mr, 
since this is still in the process of merging. Panel (b): The halo 
mass growth in dependence of the scale parameter a, fitted with 
the model given by Equation [7] to assess how typical the halo is. 




I [Gyr] 



4 INFLUENCE OF DIFFERENT ASSEMBLY 
HISTORIES 

Owing to the persistent problem of forming realistic disk 
galaxies in cosmological simulations over the past two 
decades, it became custom to choose a halo with what 
is thought to be favorable conditions for the formation of 
a disk. This comprises particularly a high spin parame- 
ter, resulting in disks with large scale lengths, and no late 
large mergers which could destroy the disk. The goal is 
to succeed in simulating at least one good disk as a basis 



Figure 8. The angular momentum evolution (top panel) and 
mass growth (middle panel) for the combined stellar and gaseous 
components, and the star formation history (bottom panel) of 
the galaxies in halos with different assembly histories. The lines 
(black, red, green and blue) correspond to the halos as shown 
in Figure |7(a)| from left to right, with the black line being our 
standard halo. 
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for further work. Howev er, it has been shown recently by 
IScannapieco et ahl l|2009l ) that disks can form in a whole 
range of halos with high or low spin parameters leading 
these authors to question the previously assumed correla- 
tion between formation history and morphology of a galaxy. 
iGovernato et al.] ()2008l ) showed that an extended disk can 
form in a simulation even after larger mergers via re-growth, 
if the progenitors are gas-rich. We therefore also apply our 
standard feedback model to three other halos with a similar 
Milky Way-type mass as the standard halo, but with differ- 
ent merging histories. The only criterion in picking these ha- 
los was their mass and a distance of at least 2 Mpc h~^ to the 
next object of equal or higher mass in order to avoid direct 
encounters at z—0. In Table [l]these are the runs denoted by 
MW_mhl_mr, MW_mh2_mr and MW_mh3_mr. The assem- 
bly histories and mock observations are shown in Figure [T] 
To characterize, how typical the chosen halos are, we fit the 
mass growth curve using 



M{a) =Moe""^a= + 



(7) 



from lWechsler et al. I l|2002l ) {a is a free parameter character- 
izing a characteristic epoch of formation). For mass accretion 
histories without large mergers, this provides a good descrip- 
tion. This is the case for our standard halo, for MW_mhl_mr 
and also, with some limitation, MW_mh3_mr, as can be seen 



in Figure 7(b) 



Figure [8] shows the angular momentum, mass growth 
and star formation history of the galaxies compared to our 
standard halo. MW_mh2_mr is a quite extreme case. Two 
large baryonic progenitors merge at z«0.1, and as is obvious 
in the image in Figure 7(a) this process is not finished at 
z—0. The plot of the mass growth of the halo shows only 
a very slight increase at these late times, after a very large 
(1.26:1) merger at around z=l. The merger coincides with 
a large starburst. 

Comparing the assembly histories of MW_mhl_mr and 
MW_mh3_mr, we would have expected MW_mhl_mr to be 
more likely to have a disk, due to the large mass increase at 
z=l for MW_mh3_mr (where the mass is more than doubled; 
this is preceded by a 3.14:1 dark matter merger). However, 
as can be seen from the angular momentum evolution in Fig- 
ure [H the assembly is more chaotic for MW_mhl_mr with 
a very large loss of angular momentum at z~1.3 coincid- 
ing with a drop in star formation. After that, even though 
the disk still grows in mass, it does not in size. At z=l it al- 
ready resembles very closely its final state. It is a low angular 
momentum system with A = 0.023, compared to A = 0.05 
for MW_mh3_mr. The latter, despite its later merger, looses 
significantly less angular momentum in its merger. The an- 
gular momentum evolution and mass growth are both quite 
comparable to our standard model except for the delayed 
assembly. While in radius it is comparable to MW_mhl_mr 
at z~l, its disk can grow due to higher angular momentum 
gas still entering the system. Also, star formation sets in 
much later for this halo, peaking at instead of 3 as in 

the standard model, due to the overall later assembly (see 
Figure [8}. This helps in keeping the bulge small. Star for- 
mation is generally slightly higher for MW_mhl_mr (with 
10 Mq yr~^ from z~l to z~0.25), leading to a low final gas 
fraction of 0.015 compared to 0.084 for MW_mh3_mr, and 
overall to a much more spheroidal shape. 



We can conclude that a quiet and early assembly history 
of the dark matter halo does not necessarily lead to a disk, 
confirming the conclusions bv lScannapieco et al.l l|2009l ) . Not 
only the time of the last major merger, but rather the onset 
of star formation and the impact particularly in loss of an- 
gular momentum seems to mainly influence the disk. This 
could be rela ted to the geometry of t he merging disks as 
suggested by IScannapieco et al] (|2009l ) . The quick growth 
of the disk in our M W_mh3_mr halo is very s imilar to what 
has been reported by IGovernato et al.l (|2008l ). 



5 THE STANDARD MODEL OVER A RANGE 
OF MASSES 

To test the performance of our standard model with respect 
to the observed scaling relations of disk galaxies, we perform 
runs using halos spanning a mass range from 1.13 x 10^^ to 
1.94x 10^^ M0 in dark matter halo virial mass. For all except 
the largest mass halo, we have performed the runs with an 
effective resolution of 2048^ ("hr" runs in our nomenclature). 
The halos were selected to have no late major mergers as well 
as relatively high spin parameters. The latter did not nec- 
essarily impact the disk. The galaxies forming in the halos 
DM_hr3 and DM_hr7 are not very disk-like. While the latter 
indeed has a rather low spin parameter with A — 0.018, the 
spin parameter of the former is A = 0.034, quite close to the 
mean value of halos. All of the chosen halos have a fairly 
quiet assembly history. The gallery of mock observations for 
these runs is shown in Figured with increasing mass from 
top left to bottom right and the standard halo in high res- 
olution being the third image in the bottom row. Table [3] 
contains the characteristic parameters. 

In our sample we have two very good cases with a strong 
rotationally supported young thin disk: DM_hr2 (Mhaio ~ 
2.52 X 10" M0) and DM_hr6 (M^aio = 7.03 x lO" Mq), 
but their masses are too low to resemble the Milky Way. 
While DM_hr2 has a fairly constant star formation history 
(SFH) of about 1.75 M© yr~^, the other halos tend to have 
a series of bursts showing the self-regulation of the feedback 
model. The surface brightness profiles of the rotationally 
supported disks are very close to exponential, which can also 
be seen from their low photometric bulge-to-disk (B/D) ra- 
tios showing a bright and prominent disk. For DM_hr6, this 
is 0.53, the lowest of the set. Dynamically, the bulge tends to 
be stronger except for DM_hr2, where (B/D)dyn=0.27 only. 
Our highest mass halo does not host a disk galaxy, but a 
very luminous object. Possibly due to the lower resolution it 
was simulated with, the disk never really grows to an appre- 
ciable size. Therefore it is not expected to follow the scaling 
relationships for disk galaxies, as will be confirmed in the 
following section. This result is not surprising, since obser- 
vationally, the most luminous galaxies are indeed elliptical 
galaxies. 



5.1 Comparison with observed scaling relations 

The most important observed scaling relations with which 
one can test simulations of disk galaxies are the TuUy-Fisher 
relation, the baryonic TuUy-Fisher relation and relations 
with respect to the angular momentum. We will discuss all 
of these in the following sections. 
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Table 3. Characteristic parameters for the resulting galaxies in the standard model applied to halos of increasing masses. 

^ galaxy mass including cold gas only, in 10^" Mq; galaxy mass including all gas, in 10^° Mq; ^ in Mq yr'^; in kpc; ^ in km s'^ 

^ in IQlO Ln 




Figure 10. Compari n g the s imulated galaxies to the observed angular momentum relation (dots), well fitted by the prediction from 
iNavarro fc Steinmet j l|2000lJ) (line). The left panel shows the angular momentum-size relation, the middle panel the true angular 
momentum of gas and stars in the simulated galaxies, and the right panel the angular momentum of the rotationally supported disk of 
young stars and gas. 



5.1.1 The angular momentum relations 

For the angular momentum, we discuss two relations: 
the angular momentum-size relation (as the "exponen- 
tial disk estimator" from Section 12. 4|) and the relation of 
the rotational velocity with the actual angular momen- 
tum of the galaxy. Our methodology to measure the scal- 
ing length, rotational velocity, and angular momentum is 
descri bed in Section 2] We compare with observed data 
from iGiovanelli et al. l|l997l ). Figure [10] shows the two 
relations. The solid line is the expected relatio n (base d 
on halo properties) from iNavarro fc Steinmet3 (|2000bh . 
jcalc = 1300(Vrot/200)^ kms-^h-^kpc. The left panel shows 
the angular momentum-size relation. Since the rotation 
curves of disks are mostly flat, this translates to a mea- 
sure for the extend of the disk. For th is we calculate j using 
j = 2RdVrot l|Mo. Mao fc Whitelll99^ . The disk sizes of our 
simulated galaxies are well within the observational scatter 



of the relation. The only outlier is, as expected, DM_hr7. 
Our disks therefore do not suffer from being too compact 
and centrally concentrated, as in earlier simulations. Look- 
ing at the middle panel, which plots the actual angular mo- 
mentum content of gas and stars in the galaxy, we see that 
even though our galaxies tend to be slightly below the ex- 
pected relation, again they fall within the observed scatter. 
DMJirl, DMJir2, MW_hr and also DMJir6 have almost no 
angular momentum deficiency, while DMJirS is quite defi- 
cient due to a large amount of counterrotating gas. If we 
only look at the actual disk consisting of cold gas and young 
stars, our data scatter right around the expected relation. 
The difference between the angular momentum as calculated 
from j = 2RdVrot, the exponential disk estimator, and the 
actual angular momentum content could be equivalent to a 
systematic difference between angular momentum measured 
by observers and angular momentum derived from simula- 
tions. This prompted us to compare these values directly, as 
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Figure 9. Mock observations of runs of halos spanning a mass 
range of an order of magnitude, from 1.13 X 10^^ to 1.94x 10^^ Mq 
in dark matter halo virial mass. The top row shows runs DMJirl, 
DM_hr2, DM_hr3 and DM_hr4, the bottom row has DMJirS, 
DM_hr6, MW_hr (the standard run in high resolution) and 
DM_mr7, the only one in this set with the medium resolution 
of 1024'' effective resolution. 
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Figure 11. The true, directly measured angular momentum of 
the simulated galaxies vs. the exponential disk estimator 2R^Vrot- 
The dotted line indicates a direct correspondence. The star sym- 
bols are for the angular momentum of the whole galaxy, the tri- 
angles for the disk (gas and young stars) only. 



shown in Figure [TT] We indeed see a clear offset for the total 
angular momentum content, while the angular momentum 
of the disk (i.e. of disk gas and young stars) is directly corre- 
lated with the 'observed' values. Two conclusions arise from 
this. First, a fair comparison between observations and sim- 
ulations with respect to angular momentum is not straight- 
forward when the simulated galaxies have large bulge com- 
ponents. This is the case for those of our galaxies which 
show the largest difference, DMJirS, DM_hr4 and DM_hr5. 





DM_hr5 > 
DM_hr4 
DM_hr5 . 



log(V™ [km s"']) 

Figure 12. The I-band (top) and baryonic (bottom) TuUy-Fisher 
relations for our halos of different masses. The bl ack points in the 
top p anel are observations from the sample of ICourteau et al.l 
(|2007|). The solid line is a fit, the dashed lines the 2a observed 
scatter. The stars are for a Salpeter IMF, the diamonds for a 
Chabrier IMF for our simulated galaxies. In the bottom panel. 



the black diamonds are observational data from lMcGaughl fcoOSi ) 
with a fit shown again by the solid line. 



Our standard feedback model is not able to produce a bulge- 
less disk and only one clearly dominated by a disk, DM_hr2 
with B/D=0.27. The second conclusion is that our actual 
disks composed of young stars and cold gas have kinematic 
characteristics of real disks. 



5.1.2 The Tully-Fisher Relation 

The Tully-Fisher relation is one of the fundamental bench- 
marks for a successful model of disk formation, and one of 
the biggest challenges of simulations, since so far it remains 
difficult to reproduce the slope and the zero-point of the 
observed relation. It relates the rotational velocity with the 
total luminosity of the galaxy and is most often measured 
in the I or K-band since this is most representative of the 
total stellar mass. We compare ou r results to observations 
compiled bv lCourteau et al] |2003)- The relation is of type 
L oc V" and their fit to the logV vs. logLi plot gives a 
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slope of 0.291 ± 0.004 and a zero point of -0.835 ± 0.039 
with a scatter ainv\L ~ 0.132. The top panel of Figure [T2l 
shows our data overplotted on this relation. Our standard 
model misses the zero point of the relation by being too 
faint by about a magnitude when a Salpeter IMF is used. 
The simulated galaxies barely fall within the scatter. The 
only exception is the galaxy in halo DM_hr2 which comfort- 
ably lies within the scatter, though still below the best fit 
relation by half a magnitude. A slight improvement can be 
achieved when using a Chabrier IMF, shown by diamond 
symbols in Figure [121 but for most galaxies the discrepancy 
is still large. 

The most fundamental relation to which the Tully- 
Fisher relati on can be trace d is the baryonic TuUy-Fisher 
relation fsee lMcGaughlfioOSi )! relating rotational velocity to 
total baryonic mass. We compare to it in the bottom panel of 
Figure [121 McGaugli (2005) finds as best fit to his observa- 
tions the relation Mbar = 50Vrot. AH simulated galaxies are 
well within the observed scatter. We seem to get a slightly 
shallower relation, though we do not really have enough data 
points for a fit and the low mass halos also might be infiu- 
enced by resolution effects. With the exception of our ex- 
pected outlier, DM_mr7, we find the agreement to be quite 
good. 

From the results of the scaling relation comparisons 
we conclude that our model is able to produce galax- 
ies which are correct structurally, but still bulge domi- 
nated and not luminous enough. Our I band stellar mass- 
to- light ratios of M/L^2 (for a Chabrier IMF) are im- 
proved compared to earlier simulations for example by 
iNavarro fc Steinmet j l|2000bl ') (who had M/L«2.5), but are 
stil l too high compar e d to analytical expectations based 
on IMo, Mao fc Whitd | 1998h in combination with the ob- 
served TuUy-Fisher relation which yield M/L~l-1.5. We 
suspect this to be the main reason for the failure in fit- 
ting the TuUy-Fisher zero point, rather than an incorrect 
rotational velocity due to peaked rotation curves. Another 
effect could be that the concentration of the dark matter 
halo is too hi gh resulting in too high rotational velocities 
l|Navarro fc Stcinmcta.2000a V 



6 INCLUDING ADDITIONAL PHYSICAL 
PROCESSES 

The physical processes playing a role in the formation and 
evolution of a galaxy are much more complex than just sim- 
ple star formation and energy feedback. In this section we 
describe a number of other relevant physical effects we tested 
to study their influence on our standard model galaxy. An 
overview of the characteristic parameters of the resulting 
galaxies can be found in Table [H and Figure [131 gives a vi- 
sual impression of the results. 



6.1 Implementation 

Before discussing the results we first briefiy describe the im- 
plementation in each case. The abbreviation in brackets is 
used later when referring to the simulations. 






* 























Figure 13. Mock observations of the galaxies run with the stan- 
dard feedback model and different additional physics effects. For 
comparison the standard model itself is shown on the left. The 
top row shows the runs with UV background, type la supernovae, 
metal-dependent cooUng and the "all in" model with metal- 
dependent cooling. In the bottom row we show the runs with 
10% kinetic feedback, mass return, the "all in" model and the 
blastwave model. 




(i) A UV background (UV) 



P [g cm"'] 



Figure 14. Metal-dependent cooling clearly reduces the amount 
of hot gas in the halo (bottom panel) compared to the standard 
run (top panel). According to our definition for the size of the 
halo, gas with p > 9.29 X 10~^* g cm^"^ belongs to the halo. 



An external ultraviolet background is included in the Gad- 
get2 cooling routine available t o us fr om Volker Springel. It 
is a modified iHaardt fc Madaul (1996h spectrum with reion- 
ization at z~6 (Da ve et al.lll99S ). 

(ii) Kinetic feedback (kin) 
In order to stimulate the generation of winds by super- 
nova explosions we test a model of splitting the supernova 
energy into a thermal and a kinetic feedback component. 
The ratio is governed by a parameter fv = Ekin/EsN, which 
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Table 4. Characteristic parameters for the resulting galaxies in the standard model with different additional physical effects. The 
standard model MW_mr_encv is shown again for comparison in the top row. 

^ galaxy mass including cold gas only, in lO'^" Mq; ^ galaxy mass including all gas, in 10^" Mq; in Mq yr~^; in kpc; ^ in km s~^; 
6 in 10^° Lq 



we set to 10%. We follo w a "momentum" approach as in 
iNavarro fc Whit3 ^M), 

where the velocity of the neigh- 
boring gas particles is changed via 



Aw,- 



2UEsN{i) 



(8) 



(iii) Type la supernovae (SNI) 

Type la supernovae arise from binary systems with a to- 
tal mass of 3-16 Mq and are particularly important for a 
delayed input of feedback energy as well as for metal en- 
richment since they are a main source f or iro n. We imple- 
ment them following IScannapieco et al.l (|2006l ). assummg a 
progenitor lifetime between 0.1 and 1 Gyr. The exact time 
for energy (and metal) injection is chosen randomly from 
this range. The relative supernova rate of type la to type II 
is esti mated from observations as 0.245 (|Cappellaro et aU 
l|l999l ), for a Milky Way type galaxy (type Sbc) with 
Lb(MW) = 2.3 X 10^° Lq). Each type la supernova is as- 
sumed to have an energy of 10^^ ergs and we do not disable 
the c ooling locally for these events f following IStinson et al.l 
l2006l ). 

(iv) Mass return (met & mSN) 

In this model we return mass to neighboring gas particles in 
the same way as energy, smoothing it via the SPH Kernel. 
This should not have a big influence, but results in stars 
of different masses and therefore different amounts of re- 
turned feedback energy. This is also required to track metal 
enrichment which we do using oxygen and i ron. For type 
II supernovae we use yields from IWooslev fc Weaver! (|1995i ) 
taking into account the dependence on metallicity for the 
iron yields. To get more realistic abundances, one has to 
include type la supernovae, since they are mostly respon- 
sible for the i ron p roduction. For this we use yields from 
iRaiteri et al.l l|l996l) which are based on [Thielemann et al] 
l|l998l ) and are independent of the progenitor mass. This is 
a separate model (mSN), since the energy of type la super- 
novae is an additional factor. 



(v) Metal-dependent cooling (mete) 
The metal abundance in the gas has a strong influence on the 
strength of the cooling. At a temperature of lO^K, primor- 
dial gas will have a cooling time w hich is about 50 times 
longe r than gas with [Fe/H]=0.5 (|Sutherland fc Dopita I 
Il993h . This resuhs in a decreased amount of hot gas 
around the halo, as shown in Figure [14] in compari- 
son to the standard run. Despite its importance, most 
current s imulations do not i nclude this effect (excep 
tions are Abadi et al.l (l2003a), Scannapieco et al.l l|2006 



lOkamoto et al.1 1|2005^ and lKawata fc GibsonI l|2005l )'). Based 



on our crude model of metal enrichment described above, 
we incl ude metal-dependent cooling b y interpolating the ta- 
bles of [Sutherland fc Dopita I (|l993D which cover, besides 
the primordial case, the abundances [Fe/II]={-3.,-2.,-1.5,- 
1., -0.5, 0,0. 5} and a temperature range of 4 < log(T) < 8.5. 
For values outside these ranges, the respective extremes are 
used. 

(vi) Blastwave feedback (bl) 
In our standard model using the SPH based method of 
smoothing energy and mass over neighboring gas particles 
with a variable-sized smoothing sphere, this size is computed 
so that the enclosed mass in the sphere is constant. This is 
numerically more sensible than our first attempt of using a 
fixed (albeit physically motivated) radius of 1.37 kpc h~^. 
However, it is not very closely connected to the underly- 
ing physics it is supposed to represent. To improve on this, 
we imple mented a vers i on of the blastwave feedback sug- 
gested bv lStinson et al.l (|2006l ). In this model, the size of the 
smoothing sphere and also the time span for which cooling 
should be turned off are not free parameters, but calculated 
throu gh th e explicit blastwave solut ion based on IChevalieil 
(|l974h and lMcKee fc OstrikeJ l|l977h . The blastwave radius 
is given by 



D iril74 1^0.32 -0.16 5-0.2 

Re = 10 -Bgi tIq Fo4 PC 



(9) 
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Figure 15. Oxygen abundance in relation to iron for our stars in 
the standard model (b lack dots) compared to observational data 
from the Milky Way bv lBensbv et al.1 1|2004| ) . The green points are 
thick disk and the red points thin disk stars. Our simulation is in 
reasonable agreement with the former, though underabundant in 
iron. 



and the timescale (i.e. the time until which cooling is dis- 
abled) is 

, -, n5-92 77,0.31 0.27 5-0.64 /^riN 

tE = 10 £51 no Po4 yr (10) 

with EsN = £^51 10^^ ergs and Po4 = 10"''Pofc"\ Po and no 
are the ambient pressure and hydrogen density. We here use 
as time ts the time until the end of the snowplow phase, 
another optio n would be the time until the hot, low density 
shell survives. IStinson et al. ! (2006) f ound no sign i ficant dif- 
ference between the two. Following Stinson et al.l (|2006l ) we 
distribute feedback energy to gas particles within Re, but 
mass and metals (if turned on) within the original smoothing 
radius. 



6.2 Results 

From the visual impression of Figure [TJ] it seems that all of 
the additional physical effects except the blastwave method 
help to make the galaxy younger and brighter. This is caused 
by a further suppression of star formation early in the as- 
sembly history, leaving more gas for later star formation. 
The blastwave method also causes such a suppression, but 
in this case we additionally have a later onset of star for- 
mation and a peak reduced to 9M0yr~^ at z~1.5. Still, 
despite much fewer stars and therefore less feedback activ- 
ity, this feedback is strong enough to result in the lowest 
galaxy mass of all cases, with only 60% of the mass in the 
standard case. Photometrically the bulge is strongly reduced 
and the galaxy also has the lowest luminosity and current 
star formation rate of all models. Dynamically it is more a 
thick disk system, though it efficiently retains its angular 
momentum. Overall, the feedback in this case appears to 
be too strong, not only suppressing early, but also current 
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Figure 16. The effect of including kinetic feedback on the surface 
brightness profile for the face-on projection of the galaxy. The 
diamonds show the profile of the standard run, the crosses the 
run with kinetic feedback. The solid line is the final fit combining 
a bulge (dotted line) and a disk (dashed line) component. 



star formation too much and therefore preventing the for- 
mation of a more dominant rotationally supported disk at 
z=0. This might be offset when including other effects like 
metal-dependent cooling. Also, a better tuning of the pa- 
rameters in this model could improve the results, but that 
has to be left to future work. 

Looking now at the differences the other, more conven- 
tional individual effects can cause, we first discuss the re- 
turn of mass. As already expected, the differences here are 
small. The galaxy is systematically slightly brighter due to a 
slightly higher star formation rate which in turn is due to a 
higher cold gas fraction. The star formation history is more 
characterized by bursts, but generally the characteristics are 
comparable. 

The model including type la supernovae is also similar. 
However, this is the only model efficiently reducing the dy- 
namical bulge. The bulge-to-disk ratio B/D is 0.18 compared 
to 0.43 in the standard model. We attribute this mainly 
to a delayed peak in the star formation history, at z~2.15 
instead of z~2.6. The height of the peak is not really re- 
duced, but overall we have more gas and more younger stars 
which consequently have formed already in the disk pro- 
genitor and therefore are rotationally supported. They are 
not young enough to effectively brighten the galaxy though, 
and photometrically the bulge is more prominent. As al- 
ready mentioned above, in the model combining type la su- 
pernovae and mass return we also implemented some basic 
metal tracking. The result is shown in Figure [TJ] for oxygen 
and iron abundances in the stars. In comparison to th e ob- 
servations of Milky Way stars bv lBensbv et al.l (|2004l ). the 
agreement with thick disk stars (green points) is reasonable, 
particularly in reproducing the trend (albeit with a large 
scatter). Compared to the thin disk stars (red points) our 
results are a little underabundant in iron. Since the model is 
rather crude this agreemen t is satisfactory and also in overall 
agreement with results bv IScannapieco et al.l l|2005h . Metal 
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enrichment is necessary for including metal-dependent cool- 
ing which will be discussed below. 

The UV background shows influence on the star forma- 
tion history only at lower redshift. The rather broad star 
formation peak of the standard mode l is sharply decreased 
at zfti 2. In agreement with results bv lNavarro fc Steinmeta 
l| 19971 ) less gas cools into the disk which is now surrounded 
by a halo of hot gas. The cold gas fraction is half of that 
in the standard model. The disk is slightly smaller but also 
brighter and more defined due to increased star formation 
particularly in the outer parts. The disk has slightly less 
angular momentum than in the standard case, and it is dy- 
namically thicker at the expense of the thin rotationally sup- 
po rted component. This again follows the trend discussed 
bv lNavarro &: Steinmeta l|l997l ). who concluded that a UV 
background makes it even more difficult to form a high an- 
gular momentum disk since it predominantly reduces the 
accretion of late infalling high angular momentum gas. How- 
ever, we find that including feedback can largely offset this 
negative effect, resulting in a slightly, but not dramatically 
smaller angular momentum. 

More interesting than the UV background is the model 
with kinetic feedback since it is the only one of our physical 
models (except the blastwave model) with an almost pure 
exponential surface brightness profile for the face-on projec- 
tion as shown in Figure [16] The outer parts of the profile are 
very similar to the standard run, but the bulge component 
in the inner 5 kpc is strongly reduced. The star formation 
history is characterized by a deep drop at z ~ 2, right after 
the first peak. The feedback generated from the first stars 
is able to heat and blow out gas, permanently reducing the 
baryon fraction within the virial radius to 0.11 (compared to 
a cosmic baryon fraction of 0.167). Later, gas is not blown 
out of the halo anymore, but its accretion onto the disk and 
therefore star formation in the disk is delayed, resulting in 
a brighter disk and a reduced mass-to-light ratio from 2.18 
to 1.93, using the Chabrier IMF. The reduced star forma- 
tion particularly at early times leads to a reduction in the 
galaxy mass by 30%, which is mostly a bulge reduction. The 
reduced bulge also leads to an increase in overall momentum 
making this our only run (besides the blastwave model and 
the high resolution run) with no angular momentum defi- 
ciency at all. The dynamical decomposition still shows a 
slightly reduced rotationally supported component and an 
increase in the thick disk. 

Finally, including metal-dependent cooling basically off- 
sets the feedback effects completely. Star formation happens 
efficiently at all times, especially very early resulting in a 
large peak of « 6OM0yr~^ at z~3. The feedback model is 
not able to regulate this at all. Star formation is still hap- 
pening efficiently at z=l and a thin young disk is present, 
but overall the system is strongly bulge dominated. Other 
effects are much needed to control particularly the early star 
formation. 



7 COMBINING THE PHYSICAL MODELS 

After studying each physical effect by itself, we now per- 
form a series of runs with what we call the "all in" model, a 
combination of the UV background, kinetic feedback, type 
la supernovae and mass return. Acting in combination, they 




Figure 17. A comparison of our standard model (black) with 
the "all in" model (red) in terms of star formation history (left 
column), dynamical rotational support (middle column) and I- 
band surface brightness profile (right column). The bottom row 
shows the same comparison, but now with the "all in" model 
including metal-dependent cooling ("allmc"). 



of course will influence each other, which makes these runs 
more realistic than the individual cases we discussed be- 
fore. We flrst discuss the result for the standard halo with 
standard resolution. Visually (see Figure [13}, the resulting 
galaxy is one of the youngest with a strongly reduced bulge 
component. This is due to a very efficient suppression of 
early star formation leaving more gas for later accretion and 
leading to the highest current star formation rate of all our 
models, 5.7 Mq yr"'^. This suppression is a new effect due to 
the combination of all the models, since, though all suppress 
star formation a little at different times, none individually 
has such a large effect. This is shown in Figure [17] in the 
left panel of the top row. The middle panel shows the distri- 
bution oi jz/jc as a measure of rotational support. Clearly 
the effect of the type la supernovae can be seen, shifting 
the emphasis from a bulge to a thick disk and also slightly 
increasing the peak of the thin disk with jz/jc > 0.85. This 
thickening of the disk can also be seen from the mock obser- 
vational image of the edge on projection. Due to the quasi 
elimination of the early star formation peak, 80% of all stars 
in the final galaxy are now formed within the disk progenitor 
and therefore do not form a spheroidal configuration. But 
particularly the stars older than about 7 Gyr do not really 
remain in a tight thin disk, but their distribution thickens, 
leading to the final result. 

In the right panel of Figure 1171 we plot the surface 
brightness profile for the face-on projection of the disk and 
again we can identify the features contributed by the indi- 
vidual physical effects. The UV background is responsible 
for the increased level outside of a radius of about 18 kpc, 
with the sharp decline at the edge of the disk at « 23 kpc. 
The kinetic feedback reduces the bulge brightness in the in- 
ner region. This is not as strong here as in Figure [TS] which 
shows the model with kinetic feedback as the only effect, 
since the combination of type la supernovae, mass return 
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and UV background together result in some overall bright- 
ening of the whole galaxy. Overall, the galaxy characteristics 
are more realistic than before. Each effect plays some role 
in this, though the kinetic feedback seems to be the most 
important. 



7.1 Including metal-dependent cooling 

We now look at the "all in" run with metal dependent cool- 
ing included. Results in comparison to the standard model 
are shown in the bottom row of Figure 1171 As we have al- 
ready seen before in the run with standard feedback and 
metal-dependent cooling as an individual effect, the early 
star formation peak is strongly increased again. Since metal 
enrichment proceeds quickly in the areas where star forma- 
tion can happen, the stronger cooling of gas with higher iron 
abundance can take over quickly as well, making gas cool 
more efficiently, form stars quicker and in turn speed up the 
enrichment. However, contrary to before, now the combined 
physical effects can regulate the increased star formation. 
The peak is narrow with a sharp drop (bottom left plot of 
Figure [T7|>. The high early star formation also leads to an 
increased fraction of stars formed in clumps again, now 37% 
compared to 20% without metal dependent cooling. In com- 
bination with a lower star formation rate in the last 4 Gyr, 
this results in the strong reduction of the thin disk peak in 
the dynamical decomposition shown in the lower middle plot 
of Figure 1171 The thin disk is rather unimportant dynami- 
cally with a mass fraction of only 17% and it is more a thick 
disk-bulge system. Photometrically, we have a fairly expo- 
nential profile with the kinetic effect of reducing the bulge 
acting in combination with a reduction of brightness in the 
outer regions of the disk. The UV background effects evi- 
dently have been canceled out by the stronger cooling. From 
this we can conclude that including other relevant physical 
effects, particularly kinetic feedback, helps to overcome the 
additional problems of increased star formation due to metal 
dependent cooling. However, a first high peak of star forma- 
tion cannot be prevented, and as a result it is much more 
difficult to produce a dominant thin, rotationally supported 
disk at z=0. 



T.2 The "all in" model in high resolution 

The "all in" model is not fully converged when applied to 
our higher resolution halos. This can largely be attributed to 
the inclusion of kinetic feedback which has a rather strong 
impact on the surrounding region of the stellar particle. In 
the same way as in Figure [TTl Figure 1181 shows the effect of 
two variations of the " all in" model for the standard high res- 
olution halo, MWJir. In the bottom row of plots, the model 
is exactly the same as described in the previous section. The 
suppression of star formation is extremely efficient, creating 
a system with a rather low mass of 3.3 x 10^" M©, only 40% 
of the stellar mass of the standard model and 60% of its 
luminosity. No thin disk is formed since there is only little 
star formation after z~1.5. The galaxy has a rather low ro- 
tational velocity and is quite thick. However, in agreement 
with the impact of kinetic feedback in the standard resolu- 
tion run, the bulge component in the photometric decompo- 
sition is reduced compared to the standard run. In order to 






Figure 18. Comparison of the standard (black lines) model with 
the "all in" model (red lines) with 3% (top panel) and 10% (bot- 
tom panel) kinetic feedback for the standard halo in high resolu- 
tion. 
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Figure 19. Mock observations of the halos with different masses 
run with the "all in" model and 3% kinetic feedback. The order 
is the same as in Figure [51 except for the plot in the lower right 
panel which shows the standard halo with 10% kinetic feedback 
in the "all in" model. 



reduce the overly strong suppression of star formation and 
facilitate the formation of a young thin disk, we reduced the 
amount of kinetic feedback from 10% to 3%. As can be seen 
in the upper row of Figure 1181 this is quite successful for 
the dynamical decomposition. We find a bulge-to-disk ratio 
of 0.57 while it is around 1 in the standard run and 1.26 
in the original "all in" model. However, 3% is too little to 
drive out substantial material from the central region of the 
galaxy and the flattening of the surface brightness proflle we 
saw in the standard resolution run with 10% kinetic feedback 
disappears completely. The disk still is substantially bright- 
ened in the outer parts leading to a lower mass-to-light ratio 
(1.7 compared to 2.06) and a more prominent thin disk than 
in the standard run. The galaxy has no angular momentum 
deficiency and, due to the higher luminosity at compara- 
ble rotational velocity, the agreement with the TuUy-Fisher 
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Figure 20. Scaling relations for the different physical models as 
shown in Table l4l 



relation is improved (though still not perfect). Since this 
seems quite promising, we apply the low kinetic energy "all 
in" model to all other halos with smaller masses than the 
standard run. The resulting mock observations are shown 
in Figure [191 Generally, for the lower mass halos the model 
is not very successful, with the exception of halo DM_hr2. 
All the others are rather spheroidal with low angular mo- 
mentum and higher mass-to-light ratios than before in the 
standard model. This shows that the "all in" model even 
with lower kinetic energy is quite dependent on the halo 
characteristics and also somewhat on resolution, and more 
work is needed to investigate particularly the behavior of 
the kinetic feedback. 



8 THE INFLUENCE OF THE ADDITIONAL 
PHYSICS ON THE SCALING RELATIONS 

Figure [50] shows the different physical models overplotted 
on the observational scaling relations. Since these runs have 
all been done with our standard halo, they scatter closely 
around the same value of rotational velocity of around 
200kms"\ Most runs (with the exception of the model with 
kinetic feedback, the bleistwave model and the high resolu- 
tion model) are slightly angular momentum deficient. This 
is mainly due to the old stars in a more spheroidal config- 
uration. In all cases there is a young stellar disk which, in 
combination with the cold gas disk, fits well on the angular 
momentum-size relation. This indicates again the difficulty 
of comparing the actual angular momentum to the exponen- 
tial disk estimator. We conclude that the additional physical 
effects generally only mildly affect the angular momentum 
content, except for kinetic feedback, which retains it much 
better than the others. However, as we have shown in our 
angular momentum study in Paper 1, feedback itself is a 
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Figure 21. TuUy-Fisher and angular momentum relation for the 
different mass halos with the " all in" model with low kinetic feed- 
back. In the upper panel, the star symbols are for a Salpeter IMF 
and the diamonds for a Chabrier IMF. The star symbols in the 
lower panel are the total angular momentum for the galaxy, the 
crosses for the disk of cold gas and young stars. 



necessary ingredient to overcome the angular momentum 
problem. 

For the TuUy-Fisher (TF) relation, we find a spread 
along the relation, but no single physical effect can be tied 
to a strongly improved agreement with the observed rela- 
tion. This is consistent with the generally similar mass-to- 
light ratios, which in all our models are between 2.5 and 3 
for a Salpeter IMF and around 2 for a Chabrier IMF. While 
the agreement is good with the baryonic TF relation inde- 
pendent on the physical model used, there is some variation 
among the models with respect to the photometric TF rela- 
tion. The lower the mass-to-light ratio, the better the agree- 
ment with the TF relation. The standard resolution "all in" 
model as well as the high resolution "all in" model with low 
kinetic feedback both agree best. With M/L=1.7 the latter 
has the lowest overall mass-to-light ratio. This leads us to 
the conclusion, that in general one physical effect might in- 
duce small changes to the galaxy, but a big impact can only 
be achieved with a combination. Our "all in" model, with 
its flaws, seems to be promising in this respect. 

The results of applying the "all in" model with low ki- 
netic energy to the lower mass halos is also in agreement 
with the conclusion drawn above. Their flt to the scaling re- 
lations are shown in Figure [21] This plot indicates a connec- 
tion between angular momentum deficiency and poor agree- 
ment with the photometric TF relation. The halos disagree- 
ing with both relations also tend to have high mass-to-light 
ratios. Only in two halos, the standard halo and DM_hr2, 
does the "all in" model result in an improved (in the latter 
case actually excellent) agreement with all scaling relations. 
Both halos are the only ones close to the observed TuUy 
Fisher-relation (top panel of Figure I2ip . They are also the 
only halos with no angular momentum deficiency at all (al- 
ready in the standard model and also in the "all in" model; 
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Figure 22. The distribution of our galaxies (r ed curves) com- 
pared witli observed data (black curves) from ICourteau et al.l 
1120071 ). both normalized by their respective number of objects. 
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Figure 23. Fraction of baryons assembled into the galaxy (fbdisk) 
vs. the ratio between the specific angular momentum of the 
disk and its surrounding halo (fj). The square marked by solid 
lines denotes the region pop ulated by early simulated galaxies 
llNavarro fc Steinmet3l2000bl) . The red colored symbols show the 
standard halo with the standard model in standard (cross) and 
high (diamond) resolution. The blue colored symbols show the 
standard halo with the "all in" model for standard (star) and 
high (triangle) resolution. 



bottom panel of Figure I21|l . They have in common the qui- 
etest merging histories of all our halos. They also have rather 
stable star formation histories (halo DM_hr2 basically lacks 
an early star formation peak). If this is the prerequisite to 
form a reasonable disk galaxy, then again the question is 
raised how the number of disk galaxies in the present day 
universe can be reconciled with the high numbers of large 
mergers in a ACDM universe. 



9 DISCUSSION AND CONCLUSIONS 

We have performed a detailed study of supernova feedback 
mechanisms in a Milky Way-type halo, but also in a range 



of other halos spanning different masses and different merg- 
ing histories, using a medium and a high resolution. Out 
of a variety of methods of distributing feedback energy we 
picked the one with the best performance to be our stan- 
dard model. The important feature for the suppression of 
the formation of a large bulge was mainly to turn off cool- 
ing locally to mimic a multiphase medium. The difference 
in the results between the distribution methods themselves 
is comparatively small. The standard model is a model with 
an energy output following an exponential law, local cooling 
turnoff and a variable smoothing length for the supernova 
energy. It was then applied to halos with different masses 
and merging histories and improved with additional rele- 
vant physical effects. The first question one might ask is how 
typical are our resulting gala xies. In Figure 1221 we compare 
them to observed data (from [c"ourteau et al. 2007f ) with re- 
spect to the most interesting parameters of disk galaxies: 
the I band luminosity, the rotational velocity and the I band 
disk scale length. Our galaxies are typical in terms of their 
luminosity. The low outlier is our lowest mass halo. Their 
rotational velocities are slightly high, but well within the 
distribution. The disk scale lengths are all somewhat high. 
Even though we attempted to measure them in the same 
way as an observer would, this is not straightforward and 
there may be systematic effects, especially due to the rota- 
tion curves. While our rotation curves for the standard halo 
are mostly flat enough that the exact measuring point does 
not matter, this is not always the case for the lower mass 
halos, especially in the "all in" model. 

From our study of different physical properties we can 
conclude that a combination of type la supernovae, some 
fraction of kinetic feedback and a UV background improves 
the result significantly towards a more realistic disk. Par- 
ticularly interesting is the kinetic feedback, which efficiently 
blows out gas from the progenitor and reduces the impor- 
tance of the central bulge making the disk more exponential. 
However, this is a very harsh way of changing the gas con- 
ditions and difficult to simply extend to higher resolution. 
In our high resolution runs, we had to reduce the kinetic 
energy fraction from 10% to 3% to achieve reasonable re- 
sults, despite the resulting loss of the flattening effect in the 
surface brightness profile. Further investigation seems nec- 
essary here to improve the stability of the model. 

We also find that it is not enough to simply reduce 
the early peak in star formation. This reduction leaves a 
lot of gas for later accretion and star formation. It reduces 
the photometric importance of the bulge, but also tends to 
thicken the disk. Dynamically, the only way to decrease the 
bulge-to-disk ratio seems to be a delay of the star formation 
peak, as in our case with the effect of type la supernovae, or 
to prevent the gas to enter the disk at a later point in time. 

Including metal dependent cooling has a significant im- 
pact and actually raises the early star formation peak again, 
making it more difficult to form a disk. In combination with 
the other effects, this can be controlled, but the result is not 
quite satisfactory yet, indicating that a more complex model 
might be necessary. However, our metal enrichment is very 
basic. It seems to be in general agreement with other simula- 
tions and also with observations, but more extensive testing 
would be necessary to study enrichment processes. Still, we 
believe that including metal dependent cooling is important 
to correctly evaluate the performance of the feedback model. 
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Figur e 24. The stellar-to-halo mass relation by iMoster et al.l 
(|2009|) in comparison with our simulated high resolution halos 
of different masses. The red and blue symbols indicate the stan- 
dard, Milky Way mass halo, while the crosses are the same halo 
in standard resolution. 



We also briefly looked at a blas twave-type fe e dback . 
Our model is not as sophisticated as in IStinson et aLll|2006l ). 
since for example we did not include feedback energy output 
depending on the lifetime of the stars. In our case this feed- 
back is extremely strong and reduces overall star formation 
too much to create a Milky Way-type galaxy. However, a 
fine-tuning of parameters could lead to improvements and 
other effects might also balance this. It would be particu- 
larly interesting to look at metal-dependent cooling in this 
respect since stronger cooling at early times needs stronger 
feedback to maintain a low star formation rate. 

Our galaxies all fit well on the angular momentum-size 
relation, so they do not suffer from the earlier problem of 
too compact, centrally concentrated disks. This is primarily 
a result of the feedback, not an effect of a high resolution (see 
Paper 1). The actual angular momentum content of gas and 
stars in our simulated galaxies is mostly within the scatter 
of observed galaxies (with the exception of some of our halos 
where no real disk was formed at all), but on the lower end 
with about 60-80% of the expected relation. The only runs 
with no angular momentum deficiency are those with kinetic 
feedback and with the blastwave model. This is alleviated 
when looking only at the angular momentum of the actual 
disk consisting of young stars and cold gas. For the total 
angular momentum of the galaxy we find a systematic offset 
when comparing directly to the angular momentum calcu- 
lated, as observers do, from j = 2RdVrot, which disappears 
when looking only at the angular momentum of the disk 
made of young stars and cold gas. A comparison between 
observations and simulations is therefore only meaningful 
for clearly disk-dominated galaxies. 



9.1 The Tully-Fisher relation as a benchmark test 



The biggest challenge for a successful disk galaxy simula- 
tion remains to fit the Tully-Fisher relation. We manage to 
fit the baryonic Tully-Fisher relation quite well, as well as 
the slope for the photometric relation, but we fail to repro- 
duce the zero-point. From this we conclude that the struc- 
ture of our simulated galaxies is realistic. Including different 
physical effects scatters the galaxies along the photometric 
relation rather than improving the fit. The "all in" model, 
with a combination of several effects, does result in an im- 
provement over the standard model. For a Salpeter IMF 
we improve from an offset factor of ~2.4 to 1.9, and for a 
Chabrier IMF from «1.8 to 1.4. However, the agreement 
is still not satisfying. It is difficult to separate out the un- 
derlying reason for the problem. It can either be that the 
luminosity is too low, or that the rotational velocity is too 
high, or a combination of both. One hint might come from 
the stellar mass-to-light ratios which in our case are about 
2 (based on a Chabrier IMF) with the "all in" runs of the 
standard halo being improved to ~1.6 both in standard and 
in high resolution. These increases again to about 2 when 
metal-dependent cooling is included. Further information 
comes from Figure [23] where we foUowH avarro fc Steinmet j 
l|2000bl ) in plotting the fraction of baryons assembled into 
the galaxy vs. the ratio between the specific angular momen- 
tum of the disk and its surrounding dark matter h alo. The 
simulated galaxies bv lNavarro fc Steinmet j (|2000bh were all 
placed within the squared region indicated by the solid lines. 
While a large fraction of the baryons were assembled into 
galaxies, they had only very little angular momentum. Our 
galaxies with their much higher angular momentum content 
are much more distributed over the whole range of the plot, 
but, with respect to the baryonic content, lie still preferably 
in the upper region, above a fraction of 30% of baryons being 
assembled into the disk. Exceptions are the lower mass ha- 
los with the "all in" model. Since typical disk galaxies have 
a mass-to-light ratio of around 1-1.5, our galaxies clearly 
contain too much mass or are too dim for the stellar mass 
they have. From semi-analytical models we know that small 
baryon fractions are needed to reach an agreement between 
the stellar mass function and the halo mass function. In 
comparison with this our baryonic mass in the disk is too 
lar ge. This is confirme d by a comparison with the results 
bv lMoster etal\ (|2009l ). who derived a stellar-to-halo mass 
(SHM) relation determined by the constraint to fit the ob- 
served SDSS stellar mass function and correlation functions. 
This SHM relation is shown in Figure 1241 together with 
the results for our high resolution simulations of different 
mass halos. The simulations with the standard model are 
shown with diamonds and the ones with the "all in" model 
with triangles. The Milky Way halo results are the red and 
blue symbols, respectively, while the crosses show the corre- 
sponding result for the standard resolution halo. Our model 
produces a systematic offset towards higher stellar masses 
which is worse for the lowest mass halos. While the "all in" 
model is able to improve on that, we know from above that 
these stars form a little bit too early and therefore are too 
dim to fit the Tully-Fisher relation. As a side note it should 
be mentioned that this curve does depend on the observed 
mass of the Milky Way dark matter halo which is still under 
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Figure 25. Panel (a): TuUy-Fisher relation with luminosities computed assuming that the bulges have the same mass-to-light ratio as 
the disk. Panel (b): Same as Figure [24l but with stellar masses excluding the bulge stars. 



discussion. Recent result s from SPSS and t he RAVE survejQ 
point to a smaller mass (|Smith et al.ll2007l ) which would de- 
crease the gap b etween the curve and our simulated data 
ijXue et a l.' 2008). 

In Figurc [25] using ad-hoc assumptions, we explore ways 
on how to improve the agreement with the TuUy-Fisher rela- 
tion and the SHM relation. Looking only at our disks made 
up of mostly young stars, they have mass-to-light ratios of 
around 1.5 for a Salpeter and even 1 to 1.3 for a Chabrier 
IMF, which is in the range of observed disk galaxies. For 
the different mass halos we therefore compute new total lu- 
minosities assuming bulge luminosities based on the same 
mass-to-light ratio as the disks. This is an attempt to model 
what the galaxies would look like if the early star forma- 
tion peak would be delayed by some mechanism, possibly 
stronger stellar feedback. (We have a small delay for the 
model including type la supernova feedback, which does re- 
duce the bulge, but the effect is not big enough.) This results 
in an excellent fit of the TuUy-Fisher relation as shown in 



Figure 25(a) for the standard model (top panel) and the "all 
in" model (bottom panel). Since with this method only the 
luminosity is adjusted, but the total stellar mass is not re- 



The Radial Velocity Experiment is measures radial velocities, 
metallicities and abundance ratios for up to a million stars in the 
Milky Way. See www.rave-survey.aip.de for more information. 



duced, it would not result in an improved agreement with 
the SHM relation. If a mechanism could be devised which 
not only would delay the early star formation, but also pre- 
vent this gas from ever forming stars, and therefore a bulge, 
for example by blowing it out efficiently or keeping it hot. 



the mass would be reduced. Figure 25(b) shows that in this 



case, when only the stellar mass of the disk is taken into ac- 
count, a good agreement with the SHM relation is achieved 
for the standard model, though the shape cannot be repro- 
duced well. In the "all in" model, where the masses are lower, 
the disk alone is not massive enough to fit the relation. An 
attempt to reproduce the TuUy-Fisher relation under the 
same assumption of taking out the bulge stars completely 
was not successful. The effect of a smaller rotational veloc- 
ity due to the reduced mass is erased by the considerably 
reduced luminosity. This could be because the disk mass- 
to-light ratios are smaller, but not small enough. However, 
for the Chabrier IMF, they are close to 1 at least for some 
of the galaxies. Another possible explanation is too much 
dark matter in the halo dominating the rotation curve and 
maintaining a high rotational velocity. This high dark mat- 
ter concentration could also be a partial explanation of the 
slightly low baryon fraction we find in some of our halos. 

We conclude that even our more complex and realistic 
"all in" model is not in complete agreement with the TuUy- 
Fisher relation, though the galaxies are within the observed 
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scatter. An exception is halo DMJir2. Improvements could 
be made in the implementation of already included effects 
(like a blastwave approach) or with the addition of new ef- 
fects, for example an efficient wind model. However, from 
the discussion above, it seems that an improved fit to both 
the TuUy-Fisher and the SHM relation is possible only when 
the formation of the bulge is largely prevented and the gas 
is kept outside the galaxy permanently. None of the physical 
effects we tested was efficient enough to be responsible for 
this. Our results indicate a possible further complication in 
the high concentration of the dark matter halo, which could 
be a basic problem of ACDM. We obtain the best results 
with two halos with the most quiet merging histories, in line 
with the typical assumptions for disk galaxies but in tension 
with ACDM's hierarchical structure. 

In summary we find the question of the formation of 
realistic disk galaxies in cosmological simulations a highly 
interconnected, complex problem that is difficult to handle 
in a controlled manner. Particularly the complex feedback 
models require a recalibration when increasing the resolu- 
tion. The answer to the pending problem of forming realis- 
tic disks does not seem to be a single physical or numerical 
effect, but rather a combination of many little effects whose 
interplay has to be scrutinized carefully and systematically 
to find the solution. 
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